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ABSTRACT 

Accreting Millisecond X-ray Pulsars like IGR J00291+5934 are important because it is possible to 
test theories of pulsar formation and evolution. They give also the possibility to constrain gravitational 
wave emission theories and the equation of state of ultra dense matter. Particularly crucial to our 
understanding is the measurement of the long term spin evolution of the accreting neutron star. An 
open question is whether these accreting pulsars are spinning up during an outburst and spinning down 
in quiescence as predicted by the recycling scenario. Until now it has been very difficult to measure 
torques, due to the presence of fluctuations in the pulse phases that compromise their measurements 
with standard coherent timing techniques. By applying a new method, I am now able to measure 
a spin up during an outburst and a spin down during quiescence. I ascribe the spin up (v su — 
5.1(3) x 10~ 13 Hzs -1 ) to accretion torques and the spin down (y s< x = —3.0(8) x 10~ 15 Hzs _1 ) to 
magneto dipole torques, as those observed in radio pulsars. Both values fit in the recycling scenario 
and I infer the existence of a magnetic field for the pulsar of B ~ 2 x 10 8 G. No evidence for an 
enhanced spin down due to gravitational wave emission is found. The accretion torques are smaller 
than previously reported and there is strong evidence for an ordered process that is present in all 
outbursts that might be connected with a motion of the hot spot on the neutron star surface. 

Subject headings: stars: neutron — X-rays: stars 



1. INTRODUCTION 

The recycling scenario (Alpar et al. 1982, Radhakr- 
ishnan & Srinivasan 1982) provides an evolutionary link 
between the young slowly rotating pulsars and the old 
fast millisecond pulsars. The evolutionary phase during 
which the pulsar is spun up happens when a neutron star 
in a binary accretes gas stripped from the donor com- 
panion. The gas is channeled onto the magnetic poles 
producing X-ray pulses modulated at the rotational fre- 
quency of the neutron star. When the accreting pulsar 
starts to spin in the millisecond range it is called an Ac- 
creting Millisecond X-ray Pulsar (AMXP). 

To establish the presence of a spin up process it is 
crucial the measurement of torques. According to ac- 
cretion theory, the excess angular momentum brought 
by the accreting gas is responsible for the acceleration 
of the neutron star rotation. However, if the angular 
momentum of the gas is not sufficiently high, the neu- 
tron star can be spun down during accretion (propeller 
regime, see the seminal works of Illarionov & Sunyaev 
1975, Ghosh & Lamb 1979 and Ustyugova et al. 2006 
for a recent study). The magnitude of the spin up/down 
is correlated with the amount of accreted matter, and 
hence with the X-ray flux. In the past, several attempts 
have been made to measure a spin up/down in several 
AMXPs, and a wealth of measurements are now available 
for at least eleven out of the thirteen known AMXPs (see 
Wijnands 2004, Poutanen 2006 and di Salvo et al. 2008). 
To accomplish this, the X-ray pulse phases are measured 
and a timing model is fitted to represent the orbital mo- 
tion, the spin frequency and its first time derivative. In 
this model one makes the assumption that the rotational 
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parameters of the neutron star (spin frequency and its 
derivatives) are coincident with the observed pulse fre- 
quency and time derivatives (see for example Galloway 
et al. 2005 and Burderi et al. 2007). 

However, large deviations from this model are usually 
seen in the timing residuals. These deviations, referred to 
as X-ray timing noise, represent some unmodeled compo- 
nent in the pulse phases which has not yet a conclusive 
explanation. Recently, it has been shown that timing 
noise in at least 6 AMXPs is strongly correlated with the 
X-ray flux (Patruno et al. 2009a). The major conclusion 
was that it is the pulse phase rather than its second time 
derivative (i.e., the pulse frequency derivative) to be cor- 
related with the X-ray flux, contrary to what predicted 
by accretion theory. A similar problem for accretion the- 
ory was seen in XTE J1807-204 (Patruno et al. 2009b), 
in which the pulse frequency derivative has no correla- 
tion with the X-ray flux whereas the pulse phases are 
strongly correlated at all timescales with the X-ray flux. 
This is particularly convincing given the presence of large 
fluctuations in the X-ray flux at different timescales. 

The first claim for a detection of an accretion torque 
in an AMXP was made by Falanga et al. (2005) for the 
AMXP IGR J00291+5934 (henceforth referred to as IGR 
J00291) during the first outburst extensively observed by 
RXTE/PCA and INTEGRAL in 2004. Differently from 
many other AMXPs, the timing residuals of IGR J00291 
look quite smooth, and the presence of timing noise is not 
as dramatic as in other AMXPs (Patruno et al. 2009a). 
Therefore at a first sight, it looks quite obvious to look for 
the presence of a spin up which manifests as a parabolic 
trend in the pulse phases and can be measured according 
to standard coherent timing techniques. However, also 
the 2004 X-ray flux of IGR J00291 shows a smooth de- 
cay in time, and no sudden flux variations are observed. 
Therefore if a correlation between the X-ray flux and the 
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pulse phases is present, it is difficult to disentangle it 
from a parabolic variation due to a true spin up, which 
is expected to be uncorrelated with the X-ray flux vari- 
ations. 

Patruno et al. (2009a) studied the 2004 outburst of 
IGR J00291 and claimed that already during this out- 
burst there is a possible correlation between flux and 
pulse phases as those seen in other AMXPs. Indeed, 
these authors questioned the presence of a torque as 
strong as that detected by Falanga et al. (2005), although 
they did not quantify the magnitude of the expected 
torque in IGR J00291. 

In 2008 the AMXP IGR J00291 went in outburst again 
(Chakrabarty et al. 2008a). This outburst was quite 
anomalous with respect to the previous one observed 
in 2004, since it showed first a faint outburst with a 
peak flux of about half the value of the 2004, and then 
went into a very low flux level phase for more than 30 
days. During this low level activity phase, the source 
was not detected by RXTE and was marginally detected 
by XMM-Newton (Lewis et al. 2010). After this period, 
a new high level activity episode was recorded: the flux 
slowly rose for about 6 days, before decaying again and 
entering into quiescence on a timescale of approximately 
one week from the outburst peak. The 2008 outburst 
has shown therefore strong flux variations that might be 
correlated with the pulse phases. 

The behavior of the 2008 outburst is also very attrac- 
tive for the purpose of testing accretion theory and pul- 
sars evolution. Thanks to the long baseline of the obser- 
vations, it is possible to follow the evolution of the spin 
parameters in IGR J00291 on a timescale of 4 years. Only 
for two other AMXPs it has been possible to accomplish 
this: SAX J1808.4-3658 (Hartman et al. 2008, Hartman 
et al. 2009, Patruno et al. 2009a) and Swift J1756-2508 
(Patruno et al. 2010). The spin of these sources is con- 
sistent with very weak or no accretion torques during the 
outbursts, and with a spin down during quiescence that 
can be interpreted as a magneto dipole torque like that 
operating in radio millisecond pulsars. 

The plan of the paper goes as follows. In Section 2 
I give a detailed summary of the observations used and 
explain the methodology applied to analyze the pulse 
phases and the X-ray flux. 

In Section 3 I perform a detailed standard coherent 
analysis of the pulse phases of IGR J00291. The assump- 
tion made in this section is that accretion theory provides 
a good description of the behavior of pulse phases and 
consequently of the neutron star spin parameters, even 
if timing noise is present. This analysis is made to verify 
whether the minimal assumptions made in accretion the- 
ory are sufficient to provide a satisfactory explanation of 
the pulse phase behavior. I discuss inconsistencies in the 
results obtained with this methodology. 

In Section 4 I study the correlations between X-ray 
flux and pulse phase variations. I propose an extended 
version of the method first appeared in Patruno et al. 
(2009a) to measure the spin frequency and accretion 
torques under the assumption that X-ray flux variations 
have an effect on pulse phases. I call this method Cor- 
relation Coherent Analysis, as opposed to the Standard 
Coherent Analysis (see Sections 2.2 and 2.3). With this 
method I am able to remove the effects of the X-ray flux 
variations from the pulse phases and measure "cleaned" 



spin parameters. 

In Section 5 I discuss the implication of the measure- 
ment of the spin period of IGR J00291. I first discuss the 
behavior of the pulsar spin during quiescence (Sec. 5.1) 
and then I consider the detection of a spin-up during the 
2004 outburst (Sec. 5.2). A discussion on the origin of 
X-ray flux and pulse phase correlations is discussed in 
Section 5.3, where I suggest a possible connection of the 
pulse phase variations with a hot spot moving on the 
neutron star surface. 

The spin-up timescale of IGR J00291 and the conse- 
quences on the lack of detected sub millisecond pulsars 
and the spin distribution of accreting pulsars are dis- 
cussed in Sections 6. 

Finally, the implications of the observed spin evolu- 
tion are summarized in the framework of the recycling 
scenario (Section 7). 

2. PULSATIONS AND X-RAY LIGHT CURVE: DATA 
REDUCTION 

I use all RXTE PCA public data for the 2004 and 
2008 outbursts of IGR J00291 (Table 1). I refer to Ja- 
hoda et al. (2006) for PCA characteristics and RXTE 
PCA absolute timing. I used all available Event 125 
/us and Good- Xenon data 2 , rebinned to 1/8192 s and 
in the 5-37 absolute channel that maximizes signal-to- 
noise-ratio (S/N). The absolute channels correspond to 
w 2.5 — 16 keV. The time of arrivals are corrected to the 
solar system barycenter (TDB timescale) by using the 
best available astrometric position reported in Rupen et 
al. (2004). An optical position has also been reported 
by Torres et al. (2008) which differs by 0.25 arcsec from 
the determination of Rupen et al. (2004). If the opti- 
cal position is used, then a shift in pulse frequency and 
frequency derivative of 4 x 10 -8 Hz and 10~ 14 Hz/s is 
expected (see for example Eq. Al and A2 in Hartman 
et al. 2008). These shifts are small enough to not affect 
the results reported in the paper. 

The light curve is folded in data chunks of different 
length, between ~ 1000 and ~ 3000 seconds, keeping 
only those with S/N>3-3.3<7, giving <1 false pulse de- 
tection per source. The presence of V709 Cas in the field 
of view of IGR J00291 has no effect on the determina- 
tion of the pulse phases. Indeed, using the harmonic 
decomposition (Boynton & Deeter 1985), only the pulse 
phases at the very frequency of IGR J00291 are mea- 
sured. The pulse amplitudes are instead strongly affected 
by the presence of V709 Cas, whereas the contamination 
cancels out when calculating the ratio between the am- 
plitudes and the la error. I detect only a significant 
number of pulsations at the fundamental frequency (v) 
and then fitted the phases with a Keplerian orbit plus a 
linear and possibly a parabolic term representing v and 
v (see Section 2.2 and 2.3 for more details). 

1 constructed the X-ray light curve using the counts 
in PCA Absolute channels 5-37 (« 2.5 - 16 keV). The 
background contribution (calculated with the FTOOL 
pcabackest) is subtracted from the total counts. 

2.1. X-ray light curves for the 2004 an d 2008 outbursts 

2 Care has to be taken when combining Events 125 /is and Good- 
Xenon data, since a rigid shift of 2~ 12 s is present in the 2004 
data due to a bug in an early version of the FTOOL xcnon2fits 
(Markwardt private communication) 
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TABLE 1 

RXTE OBSERVATIONS FOR THE 2004 AND 2008 OUTBURST 



RXTE 


Year 


Start 


End 


Outburst Length 


Program IDs 




(MJD) 


(MJD) 


(days) 




2004 


53355.85 


53342.28 


13.6 


90052- 90425 


2008 


54696.77 


54691.94 


4.8 


93013- 93435 


2008 


54742.89 


54730.51 


12.38 


94408 



Fig. 1. — X-ray light curves of the two outbursts of IGR J00291. The left panel shows the 2004 data, whereas the right one refers to the 
2008 data. The y-scale is the same in the two plot, while the x-axis is not. In 2008 two faint outbursts are observed, with peak luminosities 
about half the value of 2004. The shape of 2004 and the first 2008 outburst is remarkably similar, whereas the second 2008 outburst 
shows a bell-shape light curve. In the plot are shown only points in which pulsations are detected, which correspond to minimum counts 
of approximately 8 ct/s/PCU. 

high level activity, I prefer to treat it as two separate 
outbursts since the coherent analysis will be also given 
separately (the reason of this choice is explained in Sec- 
tion 3). 

The first weak outburst has a very similar shape as 
the 2004 outburst, but a peak flux about half the 2004 
value. The outburst appeared after 3.68 yr from the 
2004 outburst (Galloway et al. 2008). Assuming that the 
spectral shape remained the same as in 2004, the average 
mass accretion rate is: M ou tb — 6 x 10~ 10 M© yr^ 1 . This 
outburst lasted for approximately 5 days, after which a 
very low level activity was recorded (Lewis et al. 2010). 

The second weak outburst appeared after 30 days and 
showed a different shape than the 2004 and the first weak 
2008 outburst. It shows a slow rise and a rapid decay 
that create interesting variations in the X-ray flux shape 
that can be tested against the pulse phase variations. 
Assuming that the spectral shape remained the same as 
in 2004, the average mass accretion rate is: M ou tb — 
5 x 10~ 10 M^yr- 1 . 

Before the 2004 outburst, other two outbursts were de- 
tected in the RXTE/ ASM, with recurrence times of 2.80 
and 3.23 yr (Remillard 2004, Galloway et al. 2008). The 
average recurrence time seems therefore to slightly in- 
crease on the long timescale. By using all these informa- 
tion, I can calculate the long term average mass accretion 
rate in IGR J00291 for two outburst/quiescence cycles: 



The X-ray light curves for the 2004 and 2008 outbursts 
are shown in Figure 1. The first outburst has an ap- 
proximately linear decay, with pulsations detected for 14 
days, and a peak luminosity of 75 ct/s/PCU in the 2.5-16 
keV energy band, corresponding to an unabsorbed flux 
of 9.5 x 10" 10 erg/s/cm 2 (Markwardt et al. 2004) for an 
Hydrogen absorption column of Njj ~ 7 x 10 21 cm -2 and 
photon index T ~ 1.7. Falanga et al. (2005) used a comp- 
tonization model and estimated a peak bolometric flux of 
~ 2.1 x 10 -9 erg/s/cm 2 , corresponding to a bolometric 
luminosity of 6.3 x 10 36 ergs _1 , assuming a distance of 
5 kpc (the distance is taken from Falanga et al. 2005). 
In this paper I fix the distance at 5 kpc in all calcu- 
lations. The corresponding mass accretion rate, given a 
neutron star mass of 1.4M© and an assumed efficiency of 
r\ = 0.15, is therefore M ~ 2 x lO -9 M0yr -1 , according 
to the expression L^oi — M r/c 2 . The mass accretion rate 
averaged over the entire 2004 outburst length is instead 
approximately M outb ~ 7 x 10 _10 M Q yr^ 1 . Of course 
uncertainties are present in this estimate, especially be- 
cause no spectral modeling has been done. Furthermore, 
there are uncertainties in the conversion of X-ray to bolo- 
metric flux and on the assumed value of r\ for the conver- 
sion of rest mass energy into radiation, so these values 
have to be taken with care and only as orders of magni- 
tude estimates. There is also a large uncertainty in the 
distance, which has been proposed to be in the range 2-6 
kpc by several authors (Torres et al. 2008, Falanga et al. 
2005, Galloway et al. 2005). 

For the 2008 outburst, the analysis is separated for 
each of the two weak outbursts observed. Although the 
2008 outburst can be considered as a single episode of 
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Both values are in excellent agreement, and point to- 
wards an accretion rate similar to that calculated for the 
AMXP SAX J1808.4-3658, which has showed six out- 
burst/quiescence cycles with an (^Al^ ~ 10~ 11 Mp 1 
(Bildsten & Chakrabarty 2001). 
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2.2. Old Method: Standard Coherent analysis 

Standard coherent methods (e.g., Taylor 1992) are 
based on folding procedures and to % 2 minimization tech- 
niques with a model describing the time evolution of the 
pulse phase 4> (t) at the barycentric reference frame. The 
pulse phases are then fitted with a Keplerian orbit and 
a spin frequency and its first time derivative (see for ex- 
ample Patruno et al. 2009b and references therein for a 
discussion of the method). I performed this fit by using 
the standard coherent timing software TEMP02 (Hobbs 
et al. 2006). The spin frequency derivative is then as- 
sociated with the accretion torque N via the simple ex- 
pression given by accretion theory: 

N = 2ttIv su (2) 

If the orbital and astrometric components are correctly 
removed from the pulse phases, one expects to see in 
the timing residuals a set of independent values which 
are normally distributed around the zero average with 
an amplitude that can be predicted by propagating the 
Poisson uncertainties due to counting statistics. 

I stress that this method uses the assumption that 
the pulse phases vary in time only because of the rota- 
tional motion of the neutron star around its spin axis and 
around the orbit, while no other effect is taken into ac- 
count. If some other process does affect the pulse phases, 
then the Standard Coherent Analysis does not return 
realistic physical parameters and their statistical uncer- 
tainties for the accreting neutron star. 

2.3. New Method: Correlation Coherent Analysis 

In this second method, one takes into account also pos- 
sible additional effects others than the rotational motion 
of the neutron star on the observed pulse phases. A 
first attempt to apply this method in AMXPs was made 
by Patruno et al. (2009a). In that paper the authors 
considered the possible influence of the X-ray flux on 
the pulse phase, suggesting that mass accretion rate M 
induced hot spot motion dominates the observed pulse 
phase variations. 

The method I present here follows the same route, and 
searches for higher/lower reference pulse frequencies than 
the one selected in Standard Coherent Analysis by using 
the correlation between pulse phases and X-ray flux. I 
make a non-trivial assumption on the nature of the corre- 
lation: there is a linear relation between the pulse phases 
and the X-ray flux. The selected pulse frequency is then 
the one that minimizes the x 2 of the linear fit between 
phase and flux, instead of the pulse frequency that mini- 
mizes the pulse phase residuals in the fit with TEMP02. 
The reason why I choose a linear relation among all the 
possible choices is because this is the simplest law, and it 
is by no means necessarily a universal law that can be ap- 
plied in all circumstances. One has to take this assump- 
tion as the "minimal hypothesis", so that it is possible 
to verify whether under the simplest circumstances it is 
possible to already obtain results which are statistically 
better than Standard Coherent Analysis. 

Differently from Patruno et al. (2009a), that used a 
constant spin frequency model, I select the best pulse 
frequency and pulse frequency derivative instead of vary- 
ing simply the pulse frequency (see also Patruno et al. 
2009b) . I scanned 3000 pulse frequency derivatives in the 



range [-10" 12 Hzs- 1 ,+10" 12 Hzs- 1 ] and 1000 pulse fre- 
quencies for each outburst. Only for the 2004 outburst 
it is possible to detect a significant spin derivative, while 
in 2008 only non-constraining upper limits were calcu- 
lated. Since no significant pulse frequency derivative is 
detected in 2008, I refit the same data with a constant 
spin frequency model. 

By applying this technique I obtain results for all the 
outbursts which show clear improvement in the \ 2 °f the 
fit when using the Correlation Coherent Analysis rather 
than Standard Coherent Analysis. It is relevant to note 
that when using the Correlation Coherent Analysis in- 
stead of Standard Coherent Analysis, the number of pa- 
rameters used in the fit (and therefore the degrees of 
freedom) are exactly the same for the two methods, so 
there is no risk to over-fit the data when using the former 
method. 

In the following Sections 3 and 4 I discuss the re- 
sults and the implications of the two methods. If one 
applies Standard Coherent Analysis, then the physical 
consequences seem contradictory and require extraordi- 
nary explanations. If one considers Correlation Coherent 
Analysis, then the results match with a high degree of ac- 
curacy the predictions of Accretion Theory. 

3. STANDARD COHERENT ANALYSIS AND ACCRETION 
THEORY 

The spin parameters for the three outbursts are shown 
in Table 2. The orbital parameters are instead shown in 
Table 3, 4 and 5. 

According to this method of analysis the spin frequency 
between the first and the second 2008 outburst decreases 
by ~ 0.7 /i Hz. This means that during the first 2008 out- 
burst, or in between the two 2008 outbursts, a strong spin 
down must have occurred. However, a spin down during 
the outburst is not a realistic hypothesis because we de- 
tect positive spin frequency derivatives, firmly excluding 
any spin down. According to this method of analysis, 
the required spin down must have occurred in between 
the two 2008 outbursts and it needs to be of the order of 
—3 x 10 _13 Hzs _1 . This value appears to be too large to 
be explained with a propeller scenario. Indeed, at MJD 
54703, during the 30 days of low level activity, an XMM- 
Newton observation was performed, with upper limits on 
the 0.5-10 keV X-ray flux of ~ 10~ 14 erg s _1 cm -1 (Lewis 
et al. 2010), which are almost 5 orders of magnitude lower 
than the peak luminosity of the first 2008 weak outburst. 
This means that the accretion level must have been ex- 
tremely low, around 10 -14 M yr _1 or less, and this is 
not compatible with the required mass that need to be 
ejected from the system (which is of the order of a few 
10~ 10 Mq yr~ r ) to explain the ~ 0.7 /iHz spin frequency 
shift. Of course one can argue that at the onset of the 
propeller a large amount of mass is still present, but the 
X-ray luminosity is suppressed by the lack of accretion on 
the neutron star surface. However, also the production 
of X-rays in the inner accretion disk must be suppressed, 
since the X-ray luminosity observed is comparable with 
the quiescent luminosity (Campana et al. 2008; Torres et 
al. 2008). 

Also a magnetic dipole induced spin down can be con- 
sidered a quite unlikely possibility since the spin down 
between the two 2008 weak outbursts requires a magnetic 
field of at least 2 x 10 9 G and there is no reason why such 
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Fig. 2. — X-ray flux vs. pulse phase correlation for the 2004 and the two 2008 outbursts. All three outbursts follow a correlation between 
flux and phase which is consistent with being the same in all three cases, within the statistical errors. The correlation has been fit with a 
linear relation between flux and pulse phase. The pulse phase shifts by almost 0.2 cycles during 2004 and by ~ 0.1 cycle during 2008, for 
both outburst episodes. The fit parameters and the \ 2 of the fit are reported in Table 6. 

large spin down is not observed between the 2004 and the 
first 2008 outbursts. Indeed, the spin frequency at the 
end of the 2004 outburst is 598.89213111(5) Hz and with 
a spin-down of —3 x 10 _13 Hzs _1 constantly operating 
in the ~ 1336 days of quiescence, the spin frequency ex- 
pected at the beginning of the first 2008 outburst should 
have been ~ 598.8920965 Hz. This is 34^ Hz off the mea- 
sured spin frequency at the beginning of the first 2008 
outburst, and hundreds of sigma away. 

A final possibility to explain these mismatching spin 
frequencies is via glitches occurring during quiescence, 
since no signature of glitches is seen in the timing anal- 
ysis of the outbursts. Glitches in AMXPs have never 
been observed and those detected in millisecond (radio) 
pulsars are a rare phenomenon. In these latter systems, 
the magnitude on the spin frequency variation during 
a glitch is orders of magnitude smaller than in normal 
young pulsars (Cognard & Backer 2004). There is only 
one detection of a glitch in a slowly rotating accreting 
pulsars for the source KS 1947+300 (Galloway et al. 
2004). However, in this accreting pulsar the glitch re- 
sulted in an acceleration of the spin rotation, while we 
need here a deceleration. Therefore, although a glitch 
cannot be completely ruled out, it appears an unlikely 
explanation for the measured mismatching frequencies. 

I ascribe this behavior of the spin frequencies to the 
presence of timing noise as already suggested in Patruno 
et al. (2009a). A further consequence of this is that phase 
connection between the two 2008 outbursts is not a cor- 
rect procedure of analysis. Indeed, although certainly 
possible, the phase connection with the standard coher- 
ent technique ignores the presence of timing noise and 
assumes that the pulse phases are well behaved. There- 
fore the spin frequency that one finds by phase connect- 
ing the two 2008 outburst in this way is just an average 
value of the spin frequencies plus torques and contam- 



ination of timing noise and as such it is a meaningless 
quantity. However, this average value does not necessar- 
ily deviate from the true value by a large amount. This 
depends on the weighted effect of several factors, like 
the outburst length, the strength of timing noise and the 
magnitude of torques during and in between outbursts, 
that might also partially compensate each other and mis- 
lead the judgment on the validity of the results. 

Another problem is related with the magnitude of the 
spin frequency derivative detected in 2004. If one does 
not take into account the effect of timing noise, then also 
this quantity, and hence the inferred accretion torques, 
will appear higher than they really are. I give here an 
example of this effect on the measured spin frequency 
derivatives of IGR J00291. The 2004 outburst of IGR 
J00291 started on December 3, and stopped on Decem- 
ber 21. Galloway et al. (2005) analyzed the data between 
December 3 to 6 and did not detect any spin frequency 
derivative, with large upper limits due to the short base- 
line of the observation. Falanga et al. (2005) and Bur- 
deri et al. (2007) analyzed the pulsations from Decem- 
ber 7 to 21 and they both detect a consistent spin fre- 
quency derivative between 8 and 12xl0 -13 Hzs~ 1 . How- 
ever, if one splits the data differently, for example one 
analyzes the data between December 3 and 10 and be- 
tween 10 and the end of the outburst, one finds v su = 
3.6(5) x lO^Hzs- 1 and v su = 24(5) x lO^Hzs- 1 . So 
if one believes accretion theory, the torque has increased 
when the flux was lower, completely contradicting the 
expectation that a smaller amount of mass should bring 
less angular momentum and produce a smaller and not a 
higher spin frequency derivative. There might be possible 
explanations for that, i.e., considering modifications to 
the accretion theory, but taking into account the presence 
of timing noise might be a more straightforward expla- 
nation. Indeed something similar was observed in XTE 
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J 1807-294 (Patruno et al. 2009b) where the short-term 
spin frequency derivatives were behaving in a way not 
predicted by accretion theory. The conclusion was that 
the observed spin frequency derivatives were affected by 
red timing noise, so they were not true spin frequency 
derivatives. The same conclusion applies here. In the 
next Section I propose a method that takes into account 
the effect of timing noise and calculates unbiased spin 
frequencies and time derivatives. 

4. CORRELATION COHERENT ANALYSIS AND 
ACCRETION THEORY 

A possibility to explain the pulse frequency discrep- 
ancy of the previous Section is via the presence of timing 
noise in the pulse phases. To take into account this, I 
have applied the Correlation Coherent Analysis (see Sec- 
tion 2.3) to the 2004 and 2008 outburst data. 

I used the orbital solution as found in Section 3, since 
the orbit is only marginally affected by the pulse phase 
noise, which operates on completely different timescales 
than the orbital modulation. Indeed the orbital period 
of ICR J00291 is ~ 2.4 hr, while the flux changes on 
timescales of days. The very weak dependence of the or- 
bital parameters on timing noise was already noticed by 
Hartman et al. (2008) and Patruno et al. (2009b), and 
is not relevant for compact binaries in which the timing 
noise operates on timescales longer than the orbital pe- 
riod. I refer to the aforementioned papers for a detailed 
discussion of the problem and further references. 

I fit a linear relation between the pulse phases <f> and 
the X-ray flux fx- 4> — A+B- fx- The coefficients A and 
B are reported in Table 6, along with the spin frequency 
and its first derivative found minimizing the \ 2 of the 
linear fit. All coefficients are consistent with being the 
same within the statistical uncertainties. All the statisti- 
cal errors are calculated for a A\ 2 = 1 • The correlations 
for the three outbursts are shown in Figure 2, where it is 
evident that all three outbursts behave in a way which is 
consistent with being the same. Of particular relevance 
is the fact that also the second 2008 outburst follows a 
similar correlation as the 2004 and the first 2008 out- 
burst, even though the shape of this X-ray light curve 
is very different from the other two. The pulse phases 
suggest a drift of approximately 0.2 cycles in 2004 and 
0.1 cycles in the two 2008 outbursts. 

There is still some unmodeled component in the fits, 
which is particularly evident in 2004, where the pulse 
phases slightly deviate towards the end of the outburst 
and affect the goodness of the fit. However, the x 2 of 
the fits obtained with this method are statistically much 
better than those obtained with standard coherent anal- 
ysis (see last column of Table 6 and 2). The uncertain- 
ties found with this method are larger than those found 
with standard coherent analysis because here I take into 
account also the effect of timing noise in the fit. The sta- 
tistical uncertainties found in this way can be considered 
a good approximation of the true uncertainties, differ- 
ently from those obtained with standard coherent timing 
techniques. 

In a recent work, Ibragimov & Poutanen (2009) demon- 
strated that pulse phase variations have an energy de- 
pendence in the AMXP SAX J1808.4-3658. Something 
similar might also happen in ICR J00291, so it is useful 
to test whether the linear correlation and the best spin 



frequency and frequency derivative have an energy de- 
pendence. I split the data in two energy bands, a soft 
band (2-7 kcV) and an hard band (8-16 keV) and re- 
peated the entire procedure outlined above for the full 
band (2-16 keV). The best spin parameters as well as 
the coefficients A and B are consistent with being the 
same for the soft, hard and full bands. However, since 
the photon statistics is degraded when splitting the data 
in sub-bands, the statistical errors in the soft and hard 
bands are much larger than in the full band, so that a 
possible energy dependence cannot be firmly excluded. 

The long term spin frequency evolution is reported in 
Figure 3. During the 2004 outburst there is a detection of 
a spin up, while in the two weak 2008 outbursts it is only 
possible to set confidence intervals for the non-detections 
which are consistent with being the same within the sta- 
tistical uncertainties. The long term spin frequency evo- 
lution requires a constant spin down during quiescence. 
A detailed discussion follows in the next Section. 

5. LONG TERM SPIN EVOLUTION 
5.1. Pulsar spin down in quiescence 

The spin-down evolution of the neutron star in ICR 
J00291 is reported in Figure 3. The pulsar spin frequency 
requires a spin down between the end of the 2004 out- 
burst and the beginning of the first 2008 outburst. Since 
it is not possible to find a significant torque in the first 
and second 2008 outbursts, the reported spin frequencies 
refer to a constant spin frequency model. The value of 
the spin down is v s d ^ -(3 ± 0.8) x 10 _15 Hzs _1 . 

It is not possible to verify whether a spin down is also 
required between the two 2008 weak outbursts, since the 
uncertainty on the spin frequencies of the first 2008 out- 
burst is too large (~ 10~ 7 Hz). 

While other effects might also contribute to the spin 
down, the rotating neutron star magnetic field is al- 
ways present and causes a continuous emission of low 
frequency radiation. By using the force-free MHD ap- 
proximation of Spitkovsky (2006), the determination of 
v s d provides an upper limit on the dipole moment: 

M < 10 26 (l + sin 2 a )- 1/2 

(j \ 1/2 ^ 3/2 

10 45 g cm 2 J Goo Hz) 

/ - • \ 1/2 

x - - "° d G cm 3 . (3) 

V3 x 10- 15 Hz s"V w 

For the extreme values of a = 90° , 0° , and considering 
all sources of uncertainty (colatitude, period and period 
derivative), the maximum dipole magnetic field at the 
poles is 

B s d = [1-5; 2.0] ± 0.3 x 10 8 G (4) 

If the pulsar in IGR J00291 switched on as a radio 
pulsar during quiescence, then its position on the P — P 
would fall exactly in the expected region occupied by 
radio millisecond pulsars (see Figure 4). 

5.1.1. The Pulsar Magnetic Field 

Beside IGR J00291, there are only two other AMXPs in 
which it was possible to constrain the magnetic field from 
the long term spin frequency behavior: SAX J1808.4- 
3658 (Hartman et al. 2008 and Hartman et al. 2009) and 
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TABLE 2 

STANDARD COHERENT ANALYSIS: SPIN PARAMETERS FOR IGR J00291+5934 



Outburst 


Spin Frequency 


Spin Frequency Derivative 


Epoch 


X 2 /dof 




[Hz] 


[lO-^Hzs- 1 ] 


[MJD] 




2004 


598.89213045(1) 


5.6(3) 


53342.27 


309.15/88 


1st 2008 


598.89213061(8) 


12.3(4) 


54692.00 


46.39/25 


2nd 2008 


598.89213046(5) 


5.7(8) 


54730.50 


26.55/10 




1000 

Time [MJD - 53500] 



200C 



Fig. 3. — Long term spin frequency evolution of the pulsar in IGR J00291. The spin frequencies are plotted with an offset i/q = 598.8921300 
Hz. The symbols and colors are the same as those in Figure 1 and 2. The red circle and blue square are the average spin frequencies as 
measured in the first and second 2008 outburst. The cyan triangle is the spin frequency at the beginning of the 2004 outburst and the 
open cyan triangle is the spin frequency at the end of the 2004 outburst, after a spin up has taken place. The spin frequency during the 
first 2008 weak outburst has decreased with respect to 2004, and requires a spin down in quiescence of —3 X 10 5 Hzs -1 . 



TABLE 3 

2004 Orbital parameters for IGR J00291+5934 



Orbital Period [s] 

Projected Semi-major Axis [lt-ms] 

Time Of Ascending Node [MJD] 



8844.080(2) 

64.995(2) 

53345.1619277(5) 



TABLE 4 

First 2008 outburst: Orbital Parameters for IGR 
J00291+5934 



Swift J1756.9-2508 (Patruno et al. 2010). For SAX J1808 
a spin down of —5.5 x 10 _16 Hzs~ 1 was found. For Swift 
J 1756 only upper limits on the spin down were given. 

The strength of the spin down in quiescence and hence 
of the magnetic field of SAX J1808 was also revised by 
Patruno et al. (2009a), who applied a similar technique 
to that reported in this paper, and found a spin down of 
-2xl0^ 15 Hzs" 1 and a magnetic filed of ~ 2-2.8xl0 8 G, 
slightly larger than the 1.4 x 10 8 G reported in Hartman 
et al. (2009). 

The magnetic fields B of the three pulsars are there- 
fore: 



Orbital Period [s] 

Projected Semi-major Axis [lt-ms] 

Time Of Ascending Node [MJD] 



8844.069(12) 

64.987(4) 

54692.041113(3) 



IGR J00291+5934: [1.5; 2.0] ± 0.3 x 10 8 G 



SAX J1808.4-3658: [2.0; 2.8] ± 0.2 x 10 8 G 



TABLE 5 

Second 2008 outburst: Orbital Parameters for IGR 
J00291+5934 



Orbital Period [s] 8844.075(6) 
Projected Semi-major Axis [lt-ms] 64.993(5) 
Time Of Ascending Node [MJD] 54730.529224(4) 



SWIFT J1756.9-2508: [0.4; 9] x 10 8 G (95% c.l.) 



The ranges [1.5; 2.0] and [2.0; 2.8] reflect the indetermi- 
nation of the colatitude of the magnetic pole, while the 
errors on each determination is calculated by propagat- 
ing the errors on the spin parameters in Eq. (3). The B 
field of IGR J00291, SAX J1808, and Swift J1756 is per- 
fectly compatible with the minimal hypothesis that only 



TABLE 6 

CORRELATION COHERENT ANALYSIS: SPIN PARAMETERS FOR ICR .100291+5934 



Outburst 


A 




B 




Spin Frequency 


Spin Frequency Derivative [10 i3 Hzs 1 ] 


Epoch 


X 2 /dof 


2004 


-0.08 


0.01 


0.0024 


0.0003 


598.89213030(3) 


5.1(3) 


53342.27 


184.23/88 


1st 2008 


-0.10 


0.04 


0.0034 


0.0013 


598.89213061(11) 


[0; 18] (95% c.l.) 


54692.00 


46/26 


2nd 2008 


-0.08 


0.01 


0.0032 


0.0004 


598.89213070(2) 


[1; 11] (95% c.l.) 


54730.50 


19.4/11 




SsrafPjlTSS 



0.1 1 
Period (s) 



Fig. 4. — P-P diagram for radio pulsars. The cross represents 
the position of the pulsar in ICR, J00291 with period and period 
derivative as determined in Section 5.1 under the assumption that 
the observed spin down is caused by magneto dipole torques. I 
plot also the position of SAX J1808. 4-3658 (red star, Patruno et 
al. 2009a) and Swift J1756.9-2508 (blue circle, Patruno et al. 2010). 
The latter has no detected spin down during quiescence, and there- 
fore the position is determined by using an upper limit on the spin 

magneto dipole torques are at work. Under this assump- 
tion, the derived magnetic fields are exact values and not 
upper limits. There is no clear evidence for an alterna- 
tive/additional mechanism to explain the observed spin 
down other than magneto dipole torques. 

5.2. Pulsar spin up during outbursts 

According to accretion theory, a pulsar accreting from 
a disk will experience a positive torque (spin-up) when 
the magnetospheric radius r m is smaller than the coro- 
tation radius r co . The latter is defined as the position at 
which the gas in the accretion disk has the same angular 
velocity of the neutron star: 

2/3/ M xl/3 
■ 17 km ■ ( ! | t-ttt- I (5) 



lms J \\AM Qi 

For a pulsar spinning at ~ 600 Hz like IGR J00291, the 
corotation radius is at ~ 24 km, which is within approx- 
imately one stellar radii from the neutron star surface. 

The magnetospheric radius is related to the neutron 
star magnetic dipole moment [i, to the neutron star mass 
M and to the mass accretion rate M via the equation: 

/ u \ 4 / 7 
r m ~35km£^- 



,10 26 Gcm 3 



10 



-10 



M Q yr 



-1\ 2 /7 



M 



1.4 M, 



M 



1/7 



(6) 



The accretion disk mo del- dependent factor £ lies in 
the range ~ 0.1 — 1 (see e.g., Psaltis & Chakrabarty 
1999, Ghosh & Lamb 1979 and Bildsten et al. 1997). 
By assuming M = 1.4 M© and using the 2004 average 
mass accretion rate as calculated in Section 2.1 I obtain 
r m = 2 — 20 km, for 0.1 < £ < 1 . Values of r m smaller 
than 8-10 km are non-physical because of the presence 
of the hard surface of the neutron star. Nonetheless, 
given the large uncertainties in the determination of M 
(distance, accretion efficiency and bolometric flux) it is 
over-simplistic to favor values of £ closer to one and the 
calculations have to be interpreted just as order of mag- 
nitude estimates of the quantities involved. 

The condition r m = r co to enter the propeller regime 3 
is met when M ~ 10 _10 Mq yr _1 , which corresponds to 
an RXTE count rate of ~ 5 ct/s/PCU, in excellent agree- 
ment with the minimum count rate at which significant 
pulsations are detected: ~ 8 ct/s/PCU. Although this 
value depends on several additional sources of uncer- 
tainty , it shows that the magnetospheric radius really 
lies within r co during the whole period in which pulsa- 
tions are observed from 2004 to 2008. 

The maximum observed X-ray flux, and hence the 
maximum mass accretion rate as determined in Sec- 
tion 2.1 is M = 2 x 10~ 9 M© 2/r -1 . The magnetospheric 
radius at this flux is r m ~ 1.5 — 15 km, for 0.1 < £ < 1. 

To check the consistency of our results I can compare 
the maximum expected dipole magnetic filed calculated 
with accretion theory with the magneto dipole torque 
inferred independently from the spin down (Section 5.1 
and 5.1.1). Under the minimal assumption that r m — 



fi max = 2.0 x 10 23 Gcm" 
/ M 



1/2 



10" 1U M Q yr~ 



7/4 



1.4 M, 



M 



1/4 



(7) 



The value of Umax, and hence of the magnetic dipole 
field depends on the choice of £ and on the average 
outburst mass accretion rate. By choosing the least 

constraining values for £ = 0.1 and (M) = 6 x 



10~ 10 A/q yr _1 , I obtain a maximum dipole magnetic 
field at the magnetic poles of B max ~ 1.5 x 10 10 G. For 
£ = 1 the maximum B field becomes B max ~ 2.5 x 10 8 G. 

3 More recent works (like Rappaport et al. 2004 and Spruit & 
Taam 1993) show that the onset of the propeller is more compli- 
cated than simply r m = r co - However, this does not substantially 
affect the argument used in this paper. 
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These values are within the expected range for accret- 
ing millisecond pulsars and they match the spin down 
dipole magnetic field B s d as well as the upper limits of 
B m ax < 3 x 10 8 G proposed by Torres et al. (2008). 

To further check the self consistency of the accretion 
theory, it is useful to obtain a minimum magnetic field. I 
use the same argument as used in Psaltis & Chakrabarty 
(1999) and Miller ct al. (1998): pulsations are seen at 
the outburst peak, therefore the magnetic field must be 
strong enough to channel the accretion flow and enforce 
corotation of the accreted gas when the accretion rate is 
maximum. By using the least constraining mass accre- 
tion rate M = 2 x lO^M©*//- -1 (see Section 2.1) and 
using Eq. (11) in Psaltis & Chakrabarty (1999) I obtain 
a minimum magnetic field at the poles B m i n = 6 x 10 7 
G. In this calculation I have assumed the conservative 
quantities R = 10km, M = 1.4 Mq and £ = 1 (here £ is 
what is called jb in Psaltis & Chakrabarty 1999). 

Both maximum and minimum magnetic fields as in- 
ferred from accretion theory are therefore consistent with 
the spin-down magnetic dipole field obtained in Sec- 
tion 5.1. 

5.3. Motion of the Hot Spot 

In Section 4 I showed that in both 2004 and in the two 
weak outbursts of 2008 the coefficients of the linear corre- 
lation are consistent with being the same, which suggests 
some ordered process acting during each outbursts that 
provides an identical response of the phase to a fixed 
perturbation of the X-ray flux. 

Lamb et al. (2008a) suggested a moving hot spot on 
the neutron star surface as the origin of timing noise in 
AMXPs. Moving hot spots were also observed in MHD 
simulations of Romanova et al. (2003) and Romanova 
et al. (2004) although the number of simulated neutron 
star rotations is still too small to reach a firm conclusion. 

Applying the correlation coherent technique, I find 
that the pulse phases drift by approximately 0.2 cycles 
when the X-ray flux varies by a factor ~ 10 from the 
maximum to the minimum in 2004, and by 0.1 cycles in 
the two 2008 outbursts (see Figure 2). Such variations 
might be induced by a motion of the hot spot. A move- 
ment along the latitude of the neutron star will produce 
only minimal changes in phase, so a drift in longitude is 
also required. 

Under the assumption that the accretion rate tracks 
the X-ray luminosity, the magnetospheric radius r m will 
move away from the neutron star surface and approach 
the corotation radius r co as the source luminosity drops 
from the peak of the outburst down to quiescence. Dur- 
ing this process, the hot spot can move about the mag- 
netic pole depending on the colatitude of the magnetic 
axis, as it was observed in MHD simulations by Ro- 
manova et al. (2003). 

In a recent work, Bachctti et al. (2010) performed de- 
tailed MHD simulations of accreting neutron stars, and 
found that the hot spot does indeed move about the 
magnetic pole, with small colatitudes favoring a more 
pronounced motion. Although these authors could simu- 
late only a few rotational cycles of the neutron star, due 
to computational power limitations, their results clearly 
show that the accretion flow in the pulsar magnetosphere 
is a highly dynamical process and questions the tradi- 
tional picture of a static hot spot. 



A detailed discussion on the modeling of the hot spot 
movement is beyond the scope of this paper. A model 
taking into account variations of the X-ray flux and the 
hot spot motion will be presented in an accompanying 
paper. 

6. SPIN EQUILIBRIUM AND SPIN UP TIMESCALE 

A self-consistent measurement of the spin frequency 
and its long term evolution has been given, and 
by knowing the approximate duty cycle for the out- 
bursts/quiescence cycle of IGR J00291, it is possible to 
determine a timescale for the spin up of the pulsar. 

I call A^ 1 the inverse of the duty cycle for the out- 
bursts/quiescence episodes, whose value is approximately 
1% (2004 outburst length- 13 d, 2004-2008 quiescence 
length— 1350d) which is a term necessary because the 
accretion torques act only during the outbursts, while 
they are not effective during quiescence. In this first ap- 
proximate calculation I do not include the effect of the 
spin down due to magnetic dipole torques in order to ob- 
tain the shortest timescale possible. I use v s ~ 600 Hz 
and v su ~ 5 x 10~ 13 Hzs -1 , with a duty cycle A ~ 0.01. 
Therefore the spin up timescale is 

tspin-up = 4-^-A _1 = 4 Gyr (8) 

which means that the pulsar will not change its spin 
frequency significantly for a long timescale, or in other 
words that it is close to the spin equilibrium (see for 
example Campana et al. 1998). 

Furthermore, there is a further slow down on this 
timescale deriving from the spin down observed in qui- 
escence. The timescale for the spin-down (v s( i ~ —3 x 
lO-^Hzs- 1 ) is: 

ts P in-down = tt^t(1 - A) -1 ~ 8 Gyr (9) 

\Vsd\ 

which is of the same order of magnitude as the spin- up 
timescale. By combining these two timescales, under the 
hypothesis that they are both representative of the long 
term behavior of the pulsar, the long-term timescale for 
the spin evolution is: 

tspin = (-. t, £r— x) 7 Gyr (10) 

\v sd x (1- A) + u su x A J 

This timescale strengthens the suggestion that the pulsar 
in IGR J00291 is indeed close to spin equilibrium. 

6.1. The lack of Sub- Millisecond Pulsars 

A very important problem in neutron star physics is 
how to justify the absence of pulsars with spin period in 
the sub-millisecond rage. Although only a rather limited 
number of sources were known at the time, Chakrabarty 
et al. (2003) and Chakrabarty (2005) studied the spin 
distribution of AMXPs and nuclear powered pulsars and 
discovered that this was consistent with a flat distri- 
bution truncated at approximately 730 Hz with a 95% 
confidence level. Although strong observational biases 
exist for the detection of a radio pulsar above this fre- 
quency, X-ray observations taken with observatories like 
RXTE/PCA should not be affected by a significant loss 
of sensitivity at least up to 2 kHz (see Chakrabarty 2008 
and references therein). Hence it remains unexplained 
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TABLE 7 

Accretion and Nuclear Powered Millisecond Pulsars 



Nuclear-powered Pulsars 
Accretion-powered Pulsars 



Sourc6 Name 


Spin Frequency [Hz] 
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MXB 1659-298 
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4U 1636-536 


581 


IGR J00291-5934 


599 


SAX J1750.8-2900 


601 


4U 1608-52 


620 



why the maximum spin frequency known for accreting 
neutron stars is 619 Hz (Hartman et al. 2003). Hessels 
et al. (2006) discovered the fastest known radio millisec- 
ond pulsar, that spins at 716 Hz, surprisingly close to the 
cut-off spin limit of 730 Hz. If I repeat the calculation 
of the spin distribution cutoff as done by Chakrabarty et 
al. (2003), with a sample size which has doubled in the 
meanwhile (from 11 to 23 known accreting neutron star 
spins, see Table 7) I still find a cutoff of 730 Hz, but with 
a higher confidence level of 99%. 

IGR J00291 is the third fastest known accreting neu- 
tron star and among the fastest neutron stars known. It 
is therefore at the upper end of the spin distribution of 
accreting neutron stars, which is plotted in Figure 5 (see 
Watts et al. 2008,Markwardt et al. 2007 and Galloway 
et al. 2010 for references). If the spin evolution of this 
source is representative of the behavior of AMXPs, then 
it explains the existence of a cutoff of 730 Hz in the spin 
distribution of accreting pulsars. 

Several mechanisms have been proposed to explain the 
lack of sub-millisecond pulsars, but what appears cer- 
tainly true is that these spin frequencies are not limited 
by the break-up frequency of neutron stars, which sets 
in when the centrifugal force exceeds the gravitational 
pull. The break-up frequency depends on the equation 
of state of ultra-dense matter, and its determination is 
therefore of fundamental importance to understand the 
ground state of matter (Weber 2005). Basically all equa- 
tion of states allow a much higher break-up frequency 
than 730 Hz. Alternative models to explain the lack of 
sub-millisecond accreting neutron stars have been pro- 
posed, including loss of angular momentum due to emis- 
sion of gravitational radiation and the existence of mag- 
netic spin equilibrium. 

In the former case, a train of gravitational waves is 
emitted as soon as the neutron star develops a significant 
quadrupole moment. The angular momentum brought 
away by the gravitational waves slows down the pul- 
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Fig. 5. — Spin frequency distribution. The red bars represent 
10 nuclear powered pulsars whose spin frequency is measured via 
burst oscillations during thermonuclear bursts. The 13 accretion 
powered pulsars are identified by blue bars, and no accreting pulsar 
with spin below 100 Hz is counted in the sample, according to the 
definition of millisecond pulsar. Note that the red and the blue 
histograms are not overlapped, but they are added. If a pulsar is 
both accretion powered and nuclear powered, I have counted it as 
an accretion powered. The total number of neutron stars in the 
sample (accretion plus nuclear powered) is 23. 

sar rotational period, thus preventing accreting neutron 
stars to reach sub-millisecond periods. The gravitational 
wave torque is proportional to , so its effect greatly in- 
creases at high rotational frequencies. This can produce 
the sharp cutoff in the spin distribution (see Chakrabarty 
2008 for further discussions of the problem). A loss of 
angular momentum via gravitational radiation is still an 
open possibility, but the good agreement between the 
magnetic field determination via accretion torques and 
via magneto dipole spin down seems to suggest a sim- 
pler (and in this sense more likely) explanation for the 
existence of a 730 Hz cutoff. 

If the timescale for the spin-up in IGR J00291 is of 
the order of 7 Gyr (see Section 5.2), then a plausible ex- 
planation for the 730 Hz cutoff is that accreting neutron 
stars reach the spin equilibrium earlier than it is required 
to reach sub-millisecond periods. In this sense, the rea- 
son why there are no observed sub-ms pulsars might be 
a simple consequence of binary and magnetic field evolu- 
tion (see also Lamb & Yu 2005 and Lamb & Boutloukos 
2008b). 

It is possible to compare this behavior with the only 
other AMXPs with a measured spin down in quiescence: 
SAX J1808.4-3658 (Hartman et al. 2008, Hartman et al. 
2009, Patruno et al. 2009a). The overall long term spin 
frequency in SAX J1808. 4-3658 is decreasing over an ob- 
served baseline of ~ 10 years, suggesting that no sig- 
nificant accretion torque operates during the outbursts. 
Therefore also this pulsar will not significantly move in 
the spin distribution diagram on a timescale compara- 
ble with the Hubble time. Unless the neutron star was 
born with a spin already in the millisecond range, strong 
accretion torques must have spun up the pulsar in the 
past. Its current magnetic field and average mass accre- 
tion rate might be instead too small to allow a significant 
spin up during the outbursts. Therefore the spin evolu- 
tion of these two AMXPs is compatible with a scenario 
in which AMXPs evolve close to the spin equilibrium on 
a timescale shorter than it is required to spin up to the 
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sub-millisecond range. 

Before drawing a firm conclusion it is important to 
investigate the spin evolution of more accreting neu- 
tron stars, but it seems justified here to propose that, 
given the observed spin evolution of IGR J00291 and 
SAX J1808, sub-millisecond pulsars might be, at best, 
extremely rare. 

7. SUMMARY IN THE FRAMEWORK OF THE RECYCLING 
SCENARIO 

A first achievement for the recycling theory arrived 
with the discovery of the first Accreting Millisecond X- 
ray Pulsar in 1998 (Wijnands & van der Klis 1998). An- 
other fundamental step has been the recent detection of 
a millisecond radio pulsar in a position coincident with a 
a previously known quiescent neutron star X-ray binary 
(Archibald et al. 2009, Homer et al. 2006). This was a 
further evidence that accreting millisecond pulsars might 
indeed turn on as millisecond radio pulsars. 

However, there was still a missing test that needed to 
be performed before the recycling scenario could be ac- 
cepted as the correct theory of accreting neutron stars: 
the accreting pulsar is spun up, so it must be possible to 
observe accretion torques in the process of accelerating 
the neutron star rotation and spin down during quies- 
cence due to magneto dipole torques. Many claims have 
been made for a detection of an accretion torque in ac- 
creting millisecond pulsars, but none of these has been 
broadly accepted until now because of the presence of 
timing noise that affects the determination of spin and 
accretion torques when using a standard coherent timing 
analysis. 

The results presented here show that accretion torques 
are present in IGR J00291, and the spin evolution over 
4 years is entirely consistent with the prediction of the 
recycling scenario. The gas is channeled along the weak 
magnetic field lines very close to the neutron star sur- 
face, at a distance of less than 24 km from the neutron 
star center. Furthermore, a slow spin down is detected 
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when the accretion halts (or is strongly reduced), which 
I ascribe to magnetic dipole spin down as observed in 
radio pulsars. This allows the measurement of the mag- 
netic field of the neutron star, which I determine to be 
1.5 x 10 8 G < B < 2 x 10 8 ± 0.3 G, with an uncertainty on 
the two extreme values of 0.3 x 10 8 G. This value is consis- 
tent with that inferred from the accretion torques during 
the 2004 outburst. Given the large uncertainties in the 
analysis discussed in Section 5, it is still premature to 
state that the results reported in this paper finally con- 
firm the recycling scenario. However, it has been shown 
here that there is no need for new physics to explain 
the results reported. For example, there is no evidence 
for a spin down mechanism other than magneto dipole 
torques. 

There is instead strong evidence for an ordered pro- 
cess that is always present in all observed outbursts that 
might be ascribed to a motion of the hot spot on the 
neutron star surface. Finally, I find evidence for IGR 
J00291 being very close to the spin equilibrium, with the 
pulsar spin evolving on timescales of <~ 7 Gyr. These 
findings open new possibilities on the interpretation of 
a lack of detected sub-millisecond pulsars: intrinsic evo- 
lutionary processes that modify the mass transfer rate 
and the magnetic field might bring the pulsar close to 
the spin equilibrium at a much earlier stage than it is 
required to reach sub-millisecond periods. 
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